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ABSTRACT 

We present the power spectrum of solar irradiance during 1985 and 1987 
obtained from the ACRIM project from time scales of minutes to months. At 
low frequency the spectra are Lorentzian (proportional to l/(/ 2 + /„)). At 
higher frequencies they are proportional to f~%. A linear, stochastic model 
of the turbulent heat transfer between the granulation layer (modeled as a 
homogeneous thin layer with a radiative boundary condition) and the rest of 
the convection zone (modeled as a homogeneous thick layer with thermal and 
diffusion constants appropriate the lower convection zone) predicts the observed 
spectrum. 

Subject headings: sun: activity, turbulence, convection, diffusion 
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1. Introduction 

The luminosity of the sun has significant variations on time scales from minutes to 
years. Data from the NIMBUS 7 and ACRIM projects have provided us with a high quality 
time series of these variations. Some aspects of these variations can be attributed to specific 
physical processes (see Stix 1989 for an introductory review). For example, at very high 
frequencies, oscillations of the sun result in a well-understood five minute periodicity. At 
the yearly and decadal time scale there is a strong correlation between the irradiance and 
variations in the solar magnetic activity. Variations from minutes to months have a simple 
spectral form but are not well understood (Frohlich 1993). 

Frohlich (1987) has published power spectra of ACRIM data from 1980 and 1985. He 
reported that the power spectrum is flat for frequencies corresponding to time scales greater 
than one week, proportional to f~ 2 for time scales between twelve hours and one week, 
and proportional to f^ 1 at time scales down to minutes. Our studies agree that the low 
frequency spectrum is Lorentzian. We find, however, that the high-frequency spectrum is 
proportional to /~5 at time scales shorter than one day. 

Kuhn, Libbrecht, & Dicke (1988) have suggested the possibility that the interaction 
between the solar surface and deeper portions of the convection zone may cause the 
low-frequency variations. We explore that possibility in this Letter. 

If turbulent transfer in a system is dominated by eddies much smaller than the system 
size, random convective action of turbulent eddies will be analogous to the molecular 
agitation responsible for molecular diffusion (Moffatt 1983). In that approximation, 
turbulent transfer can be modeled as a stochastic diffusion process. Random forcing 
of a deterministic diffusion equation enables analytic study of the fluctuations from 
equilibrium resulting from the stochasticity of turbulent transfer in any geometry with 
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any linear boundary condition. In this paper we study the fluctuations in luminosity of 
a thin homogeneous surface granulation layer with a linear radiation boundary condition 
exchanging heat with a deep, homogeneous layer below with density and diffusion constants 
appropriate to the lower convection zone. The high and low crossover frequencies in the 
spectrum correspond to time scales of thermal and radiative equilibration of the convection 
zone, respectively. The time scales for equilibration are given by the model as a function 
of thermal and diffusion constants of the granulation layer and lower convection zone. 
Estimates of these constants with mixing length theory yields order-of-magnitude agreement 
between the crossover frequencies predicted by the model and those observed. 



2. Power Spectrum of ACRIM Data 



In Figure 1 we present the logarithm of the normalized Lomb periodograms of ACRIM 
solar irradiance data sampled during 1987 and 1985 plotted as a function of the logarithm 
of the frequency. We chose to analyze these years since they appear to represent extremes 
in the variation of solar activity at low frequencies. The low- frequency variance in the 1985 
and 1987 data are small and large, respectively. Variations in the solar irradiance at yearly 
time scales are generally agreed to be the result of variations in magnetic activity. We 
have chosen these years in order to assess the influence of magnetic activity on the power 
spectrum and distinguish its influence from that of the mechanism proposed in our model. 

Since the data were sampled at irregular intervals, simple FFT methods of estimating 
the power spectrum are only available if we average the data over some uniform time 
interval. We chose instead to use the Lomb periodogram suggested by Press et al. (1992) 
for unevenly sampled data. Above frequencies of log / = —2.0 we averaged the periodogram 
in logarithmically spaced frequency intervals of log/ = 0.01 to reduce the scatter. We 
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subtracted log S(f) by 2.5 to plot it on the same graph as the 1987 data. 

The high-frequency behavior is the same for both spectra. A f~ 2 region flattens out 
to a /~5 at frequencies greater than / Yl^Jy- ® nr observation of a scaling region at 
high frequencies disagrees with Frohlich's (1987) conclusion that the high-frequency scaling 
region is proportional to He reported this conclusion for 1985 ACRIM data, the same 
data we analyzed as part of our study. Our confidence in our interpretation lies in the 
greater resolution of our spectra, whose /~5 range is 50% larger than Frohlich's (1987). 
Large peaks appear at the orbital frequency of the satellite and its harmonics. These peaks 
are an artifact of the spectral estimation. 

The low-frequency behavior of both spectra are Lorentzian (constant at low frequencies 
and proportional to f~ 2 at higher frequencies) in agreement with Frohlich's (1987) results. 
Aside from the basic form of the spectra, there is a large variability between the crossover 
frequencies and the magnitude of the two spectra reported here and with those published 
by Frohlich (1987). This variability was also discussed by Frohlich (1987). The crossover 
frequency of the Lorentzian portion of the 1987 and 1985 data reported here are f = ■= — ^T- 

^ J 1 1 J 5 months 

and f = - — ^—rr, respectively. We interpret this variability as due to either variations in the 

•> 1 month 1 L J L J 

magnetic activity or to limitations of our model at these time scales. This Letter leaves this 
an open question. 



3. Model of Variations in Solar Luminosity 

The variations in the irradiance of the sun will be proportional to the variations in its 
surface temperature. This follows from the fact that the power emitted by the sun, modeled 
as a blackbody, F — F e = aT 4 — aT e 4 , can be well approximated by a linear dependence on 
T — T e for small departures from equilibrium. 
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Turbulent transport of heat in the convection zone of the sun can be modeled by a 
stochastic diffusion process. A stochastic diffusion process can be studied analytically by 
adding a noise term to the flux of a deterministic diffusion equation (van Kampen 1981): 

9 AT dJ 



pc- 



dt dx 



J=-a—+r ] (x,t) (2) 

where AT are the fluctuations in temperature from equilibrium and the mean and variance 
of the noise is given by 

( V (x,t)) = (3) 

( V (x, t)r](x', t')) oc a(x) (T{x)) 2 5{x - x')8(t - t') (4) 

Methods of studying transport by adding a noise term to a deterministic differential 
equation are termed Langevin methods. Langevin methods have been widely used to study 
the turbulent diffusion of a passive scalar (temperature or contaminant concentration) in 
the Earth's atmosphere and oceans (see Csanady 1980 for a review). 

A diffusion process has a frequency- dependent correlation length A = {^j-)^ (Voss & 
Clarke 1976). At very high frequencies the correlation length is smaller than the convection 
zone. In that scaling region, the fluctuations in the temperature of the granular layer of the 
sun have a power spectrum proportional to /~5 . To show this, we present an argument due 
to Voss & Clarke (1976). The Fourier transform of the driven diffusion equation presented 
above is 

AT{Kuj) = Dk^ (5) 
The frequency-dependent correlation function, c t (s,uj) =< AT(x + s,w)AT*(x,w) >, is 
then given by (Voss & Clarke 1976) 

s_ 

c(s, uj) oc ^— j- cos( j + y) (6) 

07 2 4 A 
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The power spectrum of fluctuations in the average temperature in a layer of width I will 
be a double spatial integral of the correlation function over pairs of points in the volume 
of length I and a time integral of the integrand multiplied times a factor cos(u;r) (the 
Wiener-Khintchine theorem). For low frequencies, the factor cos(wr) is one so that the 
power spectrum reduces to the spatial integral 



Since the correlation function is independent of x\ — x% for low frequencies, integration will 



At lower frequencies, the entire convection zone achieves thermal equilibrium. The 
variance in temperature of the convection zone (and therefore the irradiance) will now 
be determined by the radiation boundary condition. The fluctuating heat transport (a 
consequence of the stochasticity of turbulence) near the surface adds and subtracts heat 
from the top of the convection zone. This results in temperature and irradiance variations 
with a random walk (/~ 2 ) spectrum. 

The fluctuating input and output of heat in the f~ 2 region will cause large variations 
from equilibrium. When the temperature of the convection zone becomes larger than 
the equilibrium temperature, it will radiate, on average, more heat than at equilibrium. 
Conversely, when the temperature of the convection zone wanders lower than the equilibrium 
temperature, less heat is radiated. This negative feedback limits the variance at low 
frequencies resulting in a constant power spectrum. 

The Lorentzian portion of the spectrum has exactly the same physics as shot noise in 
an RC circuit. In shot noise, the Brownian motion of electrons gives rise to a f~ 2 spectrum 
of charge stored on the capacitor. The frequency- dependent reactance of the capacitor 
limits the variance of the charge through a transient response ^ = ^ where r D = 2nRC. 
This response is equivalent to the linear radiation boundary condition we have imposed at 




(7) 



yield the same frequency dependence. Thus, St(uj) oc lu 2. 
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the surface of the sun. In our model, two homogeneous layers with different thermal and 
diffusion constants comprise the interacting "circuit" elements. The crossover frequency is 
given by f a = 27rfi 1 c where R eq and C eq are the equivalent thermal resistance (inverse of 
the thermal conductivity) and the thermal capacitance of the parallel combination of the 
two layers, respectively. 

The model we present was solved in the context of a different problem by van Vliet, 
van der Ziel, & Schmidt (1980). They considered the temperature fluctuations in a thin 
metal film supported by a substrate. 

The geometry of the model is a thin granulation layer of width 2xl0 6 m and uniform 
density (equal to the density at the bottom of the granulation layer where most of the heat 
capacity resides) of 0.003 kg/m 3 coupled to a thick layer of uniform density representing 
the rest of the convection zone. The layers have a planar geometry for convenience in 
this simplified model. The turbulent diffusivity is estimated from mixing length theory 
as a = j^vl where v and / are the characteristic velocity and eddies sizes, respectively. 
The eddy size, I, is usually approximated as one pressure scale height. In the case of 
the granulation layer, however, the dominant eddy size is the size of the convection cell, 
approximately 2x1 6 m. The velocity at the bottom of the granulation layer is on the order 
of 1000 m/s. These estimates yield an eddy diffusivity of 10 9 m 2 /s near the solar surface. 
The thermal conductivity, a = ape, is 3xl0 7 W/m°K since c = 10 J/kg°K for a monatomic 
hydrogen gas. The velocity, thickness, and specific heat values are from a standard solar 
model presented in Stix (1992). The density values are from Bohm (1963). 

The granulation layer sits atop the rest of the convection zone. We approximate the 
density, width, and diffusivity of the remainder of the convection zone by their values near 
the bottom of the convection zone because its high density results in a concentration of heat 
capacity there and its slow diffusivity is the rate-limiting step of thermal equilibration of 
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the convection zone. The width and dominant eddy scale will both be given by 10 7 m, the 
width of the lowest pressure scale height of the convection zone. The density is estimated 
to be 0.5 kg/m 3 and the velocity is on the order of 20 m/s. These values are the arithmetic 
means of the densities and velocities at the top and bottom of the pressure scale height. 
These values yield an eddy diffusivity of 10 8 m 2 /s for the bottom of the convection zone. 
Our diffusivities agree with the estimates of Stix (1992) who quote the range of diffusivities 
in the convection zone as 10 s — 10 9 m 2 s _1 . The thermal conductivity is 3xl0 8 W/m°K. 

The equation for temperature fluctuations in space and time in the model is 
dAT(x,t) ^ d 2 AT(x,t) _ drj(x,t) 

at a[x) dx* ~ dx w 

with 

{V(x,t)) = (9) 

(r](x, t)r](x', t')) oc a{x) {T{x)) 2 5{x - x')5(t - t') (10) 

The boundary conditions are that there be no heat flow out of bottom of the convection 
zone and continuity of temperature and heat flux at the boundary separating the granulation 
layer and the deeper convection zone: 

dT 

a'—\ x=W2 = U (11) 

AT (a; = wf) = AT(x = w") (12) 
OAT ,dAT 

a -toT l *=» r-^-^-W (13) 

where W\ and u>2 are the widths of the granulation layer and deep convection zone, 
respectively and the primes denote thermal and diffusion constants of the deep convection 
zone. 

At the top of the granulation layer we impose a blackbody radiation boundary 
condition. The heat emitted from the granulation layer is dependent on the temperature 
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of the layer through the Stefan-Boltzmann law. Since solar surface temperature variations 
are small, within a linear approximation, the emitted temperature will be proportional to 
the temperature difference from equilibrium (the same approximation is commonly made in 
analytic models of climate change, i.e. Ghil 1983). At the solar surface, then, 

d AT, 



=0 = gAT(x = 0) (14) 



dx 

where g = 4:a b T^ = 2xl0 4 W/m 2o K is the thermal conductance of heat out of the sun and 
<7b is the Stefan-Boltzmann constant. 

van Vliet et al. (1980) used Green's functions to solve this model. The Green's function 
of the Laplace-transformed diffusion equation is defined by 

iujG(x, x' , iuj) — a(x) - ^—77 % ' — = 5(x — x') (15) 

ox 

where G is governed by the same boundary conditions as AT. This equation can be solved 
by separating G into two parts: G a and Gb with x < x' and x > x', respectively, where 
G a and Gb satisfy the homogeneous (no forcing) diffusion equation with a jump condition 
relating G a and G&: 

dG ai dG bl 1 

( 16 ) 



OX ox 



\x=x 



The power spectrum of the average temperature in the granualtion layer in terms of G 
is given by van Vliet et al. (1980) as: 

SATav(f) R e ( / Gi(x,x',iu)dxdx') (17) 
Jo Jo 

oc Re( / / Gib(x,x',iu)dxdx' + / / Gi a (x, x', iu)dxdx') (18) 
Jo Jo Jo Jx 

where G\ stands for the solution to the differential equation for G where the source point 
is located in the granulation layer. Two forms of G\ a and G±b are necessary for x located 
above and below x', respectively, due to the discontinuity in the derivative of G 1 created 
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by the delta function (the jump condition). The solution of G\ which satisfies the above 
differential equation and boundary conditions is 

n L (°' L ■ \,( Wl ~ x '\ ■ u( W2 \ , uf Wl ~ x \ \,( W2 \\( ■ u( x \ i a u/ x \\ 
a = ~o~K^1TL' ' — L — ' S 1 ^ U' ' — L — ' cosh( -77))( smh (^) + Y~ g cosh (^)) 

(19) 



and 



where 



G lfc = G la + -sinh(^— -) (20) 
a L 



K = smh (-i) + — cosh (-i — smh (-f) + cosh (-i) + — smh(-i)) cosh (-f) 21 
L Lg L aL' L' L Lg L L' 



and L = (-P-)s and L' = (^-)a. Performing the integration van Vliet et al. obtained 

V lb) 1 v IU1 ' ° ° 



£>AT av \f) oc Re{L (— tanh(— )(( 1) tanh(— ) — — — — + — ) 

aL' La La cosh(wi/L) L 

+(— + 7 1 - — tanh tanh -i + — — tanh -f + 1 + — tanh -± - 1 22 

a La L L g aL' L' Lg L 



For very low frequencies, 



Reducing eq. (£2] 



*nh(!£)«!*t»h(£)»£ (23) 

cosh(wi/L) - 1 _ 1 w| 
cosh(w 1 /L) 2 L 2 1 J 



1 1 

5at KI) (/) oc — — j oc (25) 
1 + J -r Jo 



which is the low-frequency Lorentzian spectrum observed in the ACRIM data. The 
crossover frequency as a function of the constants chosen for the model is 

Q a 1 . 
f n = Z ~ ~ (261 

J 2Tt(w 1 cp + w 2 c'p'{l + a f-)) 2ttw 1 w 2 c'p' 8 months v ; 

which is within an order of magnitude of the observed crossover frequencies of the 1987 and 
1985 ACRIM data, / = i — l —r- and / = = 1 ,. respectively. 

' •> 5 months •> 1 month r J 
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At low frequencies 



tanh(^) « ^, tanh(^) « 1 (27) 



L L L' 
cosh(wi/ L) — 1 1 
coshfa/L) ~ 2L2 



(28) 



then 



as observed. 

The high and low-frequency spectra meet at 



q , a . i . c'p V ,i i . cpwi .4 . . 

A = — o — s 345 -rH- 5 30 

Wipe Zfl'U'i CpO" C'p'W2 

(31) 



6 hours 

which also agrees to within an order of magnitude with the crossover frequency of the 
ACRIM data (/ = jj^). 

We have applied the same model to climate change (Pelletier 1995). The above model 
gives the fluctuations in the average temperature in a thin, homogeneous layer (Earth's 
atmosphere) coupled with another homogeneous layer with different thermal and diffusion 
constants (the ocean). As part of that work, we analyzed the Vostok ice core. We found 
the same spectral form reported here in the ACRIM data. The thermal and radiative time 
scales in that dataset were 2, 000 and 40, 000 years, respectively. Estimates of those time 
scales based upon thermal constants and diffusion constants inferred from tracer studies 
in the atmosphere and ocean matched the time scales of the crossover frequencies of the 
Vostok data well. 



4. Conclusions 
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We have presented evidence that the power spectrum of variations in solar irradiance 
exhibits four scaling regions. We presented a model originally due to van Vliet et al. (1980) 
proposed to study temperature fluctuations in a metallic film (granulation layer) supported 
by a substrate (deep convection zone) that matches the observed frequency dependence of 
the power spectrum of irradiance fluctuations. 

I wish to thank Donald Turcotte and Ed Salpeter for helpful conversations of related 
work. I am indebted to Sandy Kwan of JPL who provided me with the ACRIM data. 
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Fig. 1. — Figure 1: Logarithm of the normalized Lomb periodograms of solar irradiance in 
1987 (upper plot) and 1985 (lower plot) from the ACRIM project versus the logarithm of the 
frequency in hours -1 . The crossover frequencies for 1987 are fn = t — ^T- and f\ = , j. . 

^ •> ^ J u 5 months ■> 1 1 day 

The crossover frequencies for 1985 are f = 1 ^ onth and /i = 3-^. The 1985 spectrum is 
shifted down by log S(f) = —2.5. 



